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Abstract 

Characteristic Lie algebras of semi-discrete chains are studied. The 
attempt to adopt this notion to the classification of Darboux integrable 
chains has been undertaken. 

1 Introduction 

Investigation of the class of hyperbolic type differential equations of the form 

Uxy = f(x,y } u } u X} Uy) (1) 

has a very long history. Various approaches have been developed to look for 
particular and general solutions of these kind equations. In the literature one 
can find several definitions of integrability of the equation. According to one 
given by G. Darboux, equation JTJ is called integrable if it is reduced to a 
pair of ordinary (generally nonlinear) differential equations, or more exactly, 
if its any solution satisfies the equations of the form [I], (see also [2]) 

F(x,y,u,u x ,Uxx,...,D™u) = a(x), G(x,y,u,u y ,u yy , Dyu) = b(y), (2) 

for appropriately chosen functional parameters a(x) and b(y), where D x and 
D y are operators of differentiation with respect to x and y, u x = D x u, u xx = 
D x u x and so on. Functions F and G are called x- and ^/-integrals of the 
equation respectively. 

An effective criterion of Darboux integrability has been proposed by G. 
Darboux himself. Equation (TjQ) is integrable if and only if the Laplace se- 
quence of the linearized equation terminates at both ends. A rigorous proof 
of this statement has been found only recently [3], [1]. 
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An alternative method of investigation and classification of the Darboux 
integrable equations has been developed by A. B. Shabat based on the notion 
of characteristic Lie algebra. Let us give a brief explanation of this notion. 
Begin with the basic property of the integrals. Evidently each y-integral sat- 
isfies the condition: D y F(x, y, u, u x , u xx , D™u) = 0. Taking the derivative 
by applying the chain rule one defines a vector field X± such that 

So the vector field X\ solves the equation X\F = 0. But in general, the 
coefficients of the vector field depend on the variable u y while the solution 
F does not. This puts a severe restriction on F, actually F should satisfy 
one more equation X 2 F = 0, where X 2 = Now the commutator of 

these two operators will also annulate F. Moreover, for any X from the 
Lie algebra generated by Xi and X 2 one gets XF = 0. This Lie algebra 
is called characteristic Lie algebra of the equation ([1]) in the direction of y. 
Characteristic algebra in the x-direction is defined in a similar way. Now by 
virtue of the famous Jacobi theorem, equation flTJ is Darboux integrable if 
and only if both of its characteristic algebras are of finite dimension. In [5] 
and [B] , the characteristic Lie algebras for the systems of nonlinear hyperbolic 
equations and their applications are studied. 

In this article we will study semi-discrete chains of the following form 

ti x = f(t,t 1 ,t x ) (4) 

from the Darboux integrability point of view. Here the unknown t = t(n, x) is 
a function of two independent variables: one discrete n and one continuous x. 

It is assumed that 4f- ^ 0. Subindex means shift or derivative, for instance, 

ot x ' ' ' 

t\ = t{n + l,x) and t x = ^t(n, x). Below we use D to denote the shift 
operator and D x to denote the x-derivative: Dh(n,x) = h(n + l,x) and 
D x h(n,x) = -^h(n,x). For the iterated shifts use the subindex D^h = hj. 

Introduce now notions of the integrals for the semi-discrete chain (j3J). 
The x-integral is defined similar to the continuous case. We call a function 
F = F(x, n, t, ti, t 2 , ...) depending on a finite number of shifts x-integral of 
the chain if the following condition is valid D X F = 0. It is natural, in 
accordance with the continuous case, to call a function / = I(x, n, t, t x , t xx , ...) 
n-integral of the chain (j3J) if it is in the kernel of the difference operator: 
[D — 1)1 — 0. In other words n-integral should still unchanged under the 
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action of the shift operator DI = I, (see also [7]). One can write it in an 
enlarged form 

I(x,Tl ~\- 1, t\f f, fx , fxxi •••) I(?Gi17>)t)txi txx )•••)• (5) 

Notice that it is a functional equation, the unknown is taken at two different 
"points". This circumstance causes the main difficulty in studying discrete 
chains. Such kind problems appear when one tries to apply the symmetry 
approach to discrete equations (see [8], [9]). However the concept of the Lie 
algebra of characteristic vector fields provides an effective tool to investigate 
chains. 

Introduce vector fields in the following way. Concentrate on the main 
equation ([5]). Evidently the left hand side of it contains the variable t\ while 
the right hand side does not. Hence the total derivative of the function DI 
with respect to t\ should vanish. In other words the n-integral is in the 
kernel of the operator Y\ := D~ 1 -J=j-D. Similarly one can check that / is in 
the kernel of the operator Y 2 := D~ 2 -^-D 2 . Really, the right hand side of the 
equation D 2 I = I which immediately follows from §5§ does not depend on 
ti, therefore the derivative of the function D 2 I with respect to ti vanishes. 
Proceeding this way one can easily prove that for any natural j the operator 
Yj = D~i-^-Di solves the equation Yjl = 0. 

So far we shifted the argument n forward, shift it now backward and use 
the main equation ([5]) written as D~ l I = I. Rewrite the original equation 
(TJ1) in the form 

t-ix = g{t,t-i,t x ). (6) 
This can be done because of the condition 4f assumed above. In the 

at x 1 

enlarged form the equation D~ X I = I looks like 

I(x, 7i-l, t_i, g, g x , g xx , ...) = I(x, n, t, t x , t xx , ...). (7) 

The right of the last equation does not depend on t_i so the total derivative 
of D~ x l with respect to t_i is zero, i.e. the operator Y_i := D-^—^D^ 1 

solves the equation = 0. Moreover, the operators Y_j = D° -^—^D~i 

also satisfy similar conditions Y-jl = 0. 

Summarizing the reasonings above one can conclude that the n-integral is 
annulated by any operator from the Lie algebra L n generated by the operators 

m 

•••> Y-2, Y-i, Y-o, Yq, Yi, Y 2 , ... (8) 
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where Y = and F_o = g^~- The algebra L n consists of the operators 
from the sequence (jSJ), all possible commutators, and linear combinations 
with coefficients depending on the variables n and x. Evidently equation 
(HI) admits a nontrivial n-integral only if the dimension of the algebra L n is 
finite. But it is not clear that the finiteness of dimension L n enough for exis- 
tence of n-integrals. By this reason we introduce another Lie algebra called 
the characteristic Lie algebra of the equation (jlj). First we define in addi- 
tion to the operators Yi, Y 2 , ... differential operators Xj = -J^— for j = 1,2,.... 

The following theorem allows us to define the characteristic Lie algebra. 

Theorem 1.1. Equation ^ admits a nontrivial n-integral if and only if the 
following two conditions hold: 

1) Linear envelope of the operators {Yj}f is of finite dimension, denote this 
dimension N ; 

2) Lie algebra L n generated by the operators Yi, Y 2 , Yn, Xi, X 2 , X^ is 
of finite dimension. We call L n the characteristic Lie algebra of §4ty. 

Remark 1.2. It is easy to prove that if dimension of {Yj}f is N then the 
set {Yj}i constitute a basis in the linear envelope of {Yj}™ . 

2 Characteristic Lie Algebra L n 

Study some properties of the characteristic Lie algebra introduced in the the- 
orem 11.11 above. Begin with the proof of the remark 1 1 .21 It will immediately 
follow from the following lemma. 

Lemma 2.1. If for some integer N the operator Y/v+i is a linear combination 
of the operators with less indices: 

Y/v+i = ctiFi + a 2 Y 2 + ... + a N Y N (9) 

then for any integer j > N, we have a similar expression 

Yj = ftY x + (3 2 Y 2 + ... + (3 N Y N . (10) 

Proof. Due to the property Y^+i = D^Y^D, we have from ([9]) 

Y N+2 = D-\ ai )Y 2 + D-\a 2 )Y 3 + ... + D- l {a N ){a x Y x + ... + a N Y N ). (11) 

Now by using induction one can easily complete the proof of the lemma. 
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Lemma 2.2. The following commutativity relations take place: 



[y , r_ ] = o, [y , y 1 ] = o, [y_o, = o. 

Proof. The first of the relations is evident. In order to prove two others 
find the coordinate representation of the operators Y\ and Y_i acting in the 
class of locally smooth functions of the variables x, n,t,t x ,t xx , ... . By direct 
computations 

i d 

YJ = D~ l —DI 
dti 

= D-^IfaJJ,,...) 

(12) 

one gets 

Y 1 = — + D- 1 (^-) — + D- 1 (^) — + D' 1 ( 9fxx ) 9 + ... . (13) 
dt \dtiJ dt x \dtiJ dt xx V dtx J dt xxx 

Now notice that all of the functions /, f x , f xx ,... depend on the variables 
ti, t, t x , t xx , ... and do not depend on t 2 hence the coefficients of the vector 
field Y\ do not depend on t\ and therefore the operators Y\ and Yq commute. 
In a similar way by using the explicit coordinate representation 

1 dt (<9£_i) dt x (<%_i) dt xx (cft_i) cft xxx ^ ^ 

one can prove that [YLo, ^-i] = 0. 

The following statement turned out to be very useful for studying the 
characteristic Lie algebra L n . 

Lemma 2.3. Suppose that the vector field 

Y = a(0)d t + a(l)d tx + a(2)d txx + (15) 
where a x (0) = solves the equation [D X ,Y] = 0, then Y = a(0)d t . 
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The proof is based on the following formula 

[D x , Y] = (a x (0) - a(l))d t + (a x (l) - a(2))d tx + ... . (16) 

So if a x (0) = 0, then a(l) = 0, but if a x (l) = then a(2) = and hence 
a(j) = for all j > 0. 

In the formula ffl2l) we have already given an enlarged coordinate form of 
the operator Y\. One can check that the operator Yi is a vector field of the 
form 

Y 2 = D^QWfidt. + D-\Y x {f x ))d txx + D-\Y x {f xx ))d txxx + ... . (17) 

It immediately follows from the equation Y2 = D~ l Y\D and the coordinate 
representation (fT2l) . By induction one can prove similar formulas for arbitrary 

J- 

Y j+1 = D-\Y 3 {f))d tx + D-\Y 3 {f x ))d txx + D-\Y 3 {f xx ))d txxx + ... . (18) 

Lemma 2.4. For the operators D x , Y\, K_i considered on the space of smooth 
functions of t,t x ,t xx , ... the following commutativity relations take place: 

[D x ,Y 1 ]=pY 1 , [D x , F_i] = qY-i, 

where p = -D- l {^) and q = -D{-^). 

Proof. Recall that 

«-l + ^)K +ir, (^)£ + -- (19) 

We find [D x , Yi] by using ffToT) as 

[D x , Y x ] = -D-\f tl )d t + D-\D x (f tl ) - f xtl )d tx + ... . (20) 
For arbitrary function H, we have 

d 

[D x ,d tl ]H(t,t 1 ,t x ,t xx , ...) = D x H tl — —D X H 

Oh 

d 

= (H ttl t x + H tltl ti x + ...) — —(H t t x + H tl t lx + ...) 
= -H t J tl . ' (21) 
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By taking H = f and H = f x one gets [D x ,d tl ]f = -f t J tl , [D x ,d tl ]f x 
— fxtiftn an d so on - We insert these equations into (1201) to find 

In a similar way one can prove that [D x , YLi] = — £)^^^-^Y_i. 

Let us prove theorem 1.1. Suppose that there exists a nontrivial n- integral 
F = F(t, t x , t\jsr\) for the equation (j3J), here ty\ = D x t for any j > 0. 
Then all the vector fields from the Lie algebra M generated by {Yj,Xk} for 
j = 1,2, ... and = 1, where A^2 is chosen arbitrarily satisfying N2 > N 
annulate F . We will show that dim M < 00. Consider first the projection of 
the algebra M given by the operator Pjy: 

-1 00 -IN 

i=-N 2 i=0 i=-N 2 i=0 

Let L n (N) be the projection of M. Then evidently the equation ZqF = 
is satisfied for any Z in L n (N). Evidently, dim L n (N) < 00. Let the set 
{Z 01 , Z 02 , Z 0Nl } forms a basis in L n (N). One can represent any Z in 
L n (N) as a linear combination 

^0 = ociZ 01 + a 2 Z 02 + ... + a Nl Z 0Nl . (24) 

Suppose that the vector fields Z, Z\, Z^ x in M are connected with the 
operators Z , Z 01 , Z 0Nl in L n (N) by the formulas Pn(Z) = Z Q , P N (Zi) = 
Z i, P N (Z Nl ) = Z 0Nl . We have to prove that 

Z = aiiZt + a 2 Z 2 + ... + a Nl Z Nl . (25) 

In the proof, we use the following lemma. 

Lemma 2.5. Let Fx = D X F and F is an n-integral. Then for each Z in M 
we have ZF\ = 0. 

Proof. It is easy to check that F\ is also an n-integral, really DF\ = 
DD X F = D X DF = D X F — Fx. It was shown above that any Z in M 
annulates n-integrals. 
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Apply the operator Z — a,\Z x — a 2 Z 2 — ■■■ — ol Nx Z Nx to the function F x = 
Fi(t, t x , t xx , t[jy+i]), 



(Z - a x Z x - a 2 Z 2 - ... - a Nl Z Nl )F 1 = 0. (26) 
We can write fl26l) as 



(Z - a x Z QX - a 2 Z 02 - ... - a Nl Z mi )F x + (X(N + 1) - a x X x {N + 1) 
- a 2 X 2 (N + 1) - ... - a Nl X Nl (N + l))—?—Fi = 0, 

Ot[N+l] 

(27) 

where X(N + 1), X X (N + 1), ...,X Nl (N + 1) are the coefficients before d t[N+1] 
of the vector fields Z, Z x , Z 2 , Z Nl . The first summand in (12TI) vanishes by 
(T2U). In the second one the factor 75—^ — F x = -^f—F is not zero. So we have 

X(N + 1) = aiXi(iV + 1) + a 2 X 2 (N + 1) + ... + a Nl X Nl (N + 1). (28) 
Equation (1281 shows that 

P N+1 {Z) = a 1 P N+1 (Z 1 ) + a 2 P N+1 (Z 2 ) + ... + a Nl P N+1 (Z Nl ). (29) 



So by applying mathematical induction, one can prove the formula ( l2ol) . Thus 
the Lie algebra M is of finite dimension. Now construct the characteristic 
algebra L n by using M. Since dim M < 00, the linear envelope of the 
vector fields {^-j-! is of finite dimension. Choose a basis in this linear space 
consisting of Y x , Y 2 , Yr for K < N < N 2 . Then the algebra generated by 
Fi, 3^2, Yk-i X x , X 2 , ■■■,X k is of finite dimension, because it is a subalgebra 
of M. This algebra is just characteristic Lie algebra of the equation 

Suppose that conditions (1) and (2) of the theorem 11.11 are satisfied. So 
there exists a finite dimensional characteristic Lie algebra L n for the equation 
Oil). Show that in this case equation (j4j) admits a nontrivial n-integral. Let 
N x is the dimension of L n and iV is the dimension of the linear envelope of 
the vector fields {Yj}'jL 1 . Take the projection L n (N 2 ) of the Lie algebra L n 
defined by the operator P N2 defined by the formula (l23l) with N 2 instead of N. 

-1 N 2 

Evidently, L n (N 2 ) consists of the finite sums Z = x(i)d ti + x(i)d tli] 

i=-N i=0 

where N = N x — N 2 . Let Z ox , Z ^ 1 form a basis in L n (N 2 ). Then we have 
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Ni = N+N 2 equations Z jG = 0, j = 1, N\, for a function G depending on 
iV+A^+l = Ni+1 independent variables. Then due to the well-known Jacobi 
theorem, there exists a function G = G(t-N, t-N+i, t-x, t, t x , t xx , i[jv 2 ])> 
which satisfies the equation ZG = for any Z in L n . But really it does not 
depend on t-N, t_i because XiG = 0, X 2 G = 0, X^G = 0. Thus 
the function G is G = G(t, t x , t xx , t[jv 2 ])- Such a function is not unique 
but any other solution of these equations, depending on the same set of the 
variables, can be represented as h(G) for some function h. 

Note one more property of the algebra L n . Let 7r be a map which assigns 
to each Z in L n its conjugation D~ 1 ZD. Evidently, the map 7r acts from the 
algebra L n into its central extension L n Q){X N+ i}, because for the generators 
of L n we have D~ x YjD = Y j+1 and D^XjD = X j+1 . Evidently, [X N+1 , Yj] = 
and [Xn + i,Xj] = for any integer j < N. Moreover X^+iF = for 
the function G = G(t, t x , i[jv 2 ]) mentioned above. This fact implies that 
ZG\ = for G\ = DG and for any vector field Z in L n . Really, for any Z in 
L n one has a representation of the form D~ 1 ZD = Z + XX^ + i where Z in 
L n and A is a function. So 

ZG X = ZDG = D{D- X ZDG) = D(Z + XX N+1 )G = 0. (30) 

Therefore G\ = h(G) or DG = h(G). In other words function G = G(n) 
satisfies an ordinary difference equation of the first order. Its general solution 
can be written as G = H(n, c) where if is a function of two variables and c 
is an arbitrary constant. By solving the equation G = H(n, c) with respect 
to c one gets c = F(G, n). The function F = F(G, n) found is just n-integral 
searched. Actually, DF(G, n) = Dc = c = F(G,n). So DF = F. This 
completes the proof of the theorem 1.1. 

3 Restricted classification 

Different approaches are known to classify integrable nonlinear differential 
(pseudo-differential) equations. One of the most popular and powerful is 
based on higher symmetries. For the first time the theoretical aspects of this 
method have been formulated in the famous paper by N. Kh. Ibragimov and 
A. B. Shabat [TTj . Several classes of nonlinear models were tested by this 
method in [12], [13]. The symmetry approach allowed R. I. Yamilov to find all 
integrable chains of the Volterra type [H]: u t {n) = f(u(n— 1), u(n), u{n + \)). 
The consistency approach to the classification of integrable discrete equations 
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has been studied by Adler, Bobenko and Suris in [T5]. Classification based 
on the notion of the recursion operator is studied in [To]-[l~""]. 

In this paper we undertake an attempt to adopt the notion of characteris- 
tic Lie algebra to the problem of classification of Darboux integrable discrete 
equations of the form (j3J). The classification problem consists of describing 
of all chains admitting finite dimensional characteristic Lie algebras in both 
directions. Actually the problem of studying the algebra generated by the 
operators (|HJ) seems to be quite difficult. That is why we will start with a 
very simple case. 

Formulation of the problem. Study the problem of finding all of the 
equations (jlj) for which the Lie algebra generated by the operators Y x and 
Yli is two dimensional. 

Let us denote Y\-\ = \Y X ,Y^ X ]. We require that the following relation 
takes place Y\-\ = XY% + fiY- X . It follows from the explicit formulas ([TBI 
and (j!4p that the vector field Y\-\ does not contain a summand with the 
term hence fi = —A. The commutators of the basic vector fields with 
the operator of the total derivative admit simple expressions (see lemma 2.4 
above). Evaluate the commutator [D x , Y x _ x ] 

= \Y x ,qY_ x \-\Y- X ,pY x \ 

= Y x {q)Y_ x + qY x> _ x -Y_ x {p)Y x +pY Xt _ x 

= (p + q)Y h _ x + Y x (q)Y_ x -Y_ x (p)Y x . 

Remind that due to the reasoning above a coefficient A = X(n, x) should exist 
such that 

Y X - X = X(Y X — Y_i). (31) 

The problem is in finding of / in the equation t Xx = f(t, t x , t x ) for which the 
constraint (13*11) is valid. 

Commute each side of the equation ([ST]) with the operator D x 

[D X ,Y X ^ X ] = [D X ,XY X ]-[D X ,XY_ X ] 

= (p + q)X(Y x -Y_ x )+Y x (q)Y_ x -Y_ x (p)Y x 
= D x (\)Y x + XpY x -D x (X)Y_ x -\qY_ x . 

Compare now two different expressions for the commutator. This gives rise 
to conditions 

qX-Y_ x (p) = D x (X) , pX-Y x (q) = D x (X), 



10 



which form an over-determined system of equations for unknown A, it should 
satisfy simultaneously two equations. By solving them with respect to A and 
D X (X) we obtain equations 

A= y-.(p)-y,(») Dx[x) ^ym-vy-M (32) 

q — p p — q 

which immediately yield 

D j Y^{p)-yM \ = pY-i(p) - qYM) ^ (33) 
\ q — p J q — p 

First of all note that this equation contains both / and its inverse g. Ex- 
clude g. Recall that t lx = f(t,ti,t x ) and t x = f(t-x,t,t-i x ) where t_ lx = 
g(t,t_i,t x ). Differentiating the identity t x = f(t-x,t,g(t,t-x,t x )) with re- 
spect to t-x one gets 

D - 1 (|( Ml , ( j) +fl - 1 (|( Ml , il ))^- = (34) 

which implies that 

9t-, = -D- l {j-), (35) 

ft 

so D(g t _ 1 ) = Let us write the equation explicitly. Evaluate first 

Jt x 

Yx(q) and Y_x(p), where p = —D^ 1 (f t ) and q = 

Yx(q) = {dt + D-\f tl )d tx + D-\f xtl )d txx + ...}^- 

= (36) 

J l-I </ I'X '*' 



11 



By inserting these equations into (J33D one gets a long expression 



D 



(^(/„)) -((4) +D ->Cfi,)(4 



^(/^((^(/.j), J - £((£), + n-\M{jt 



A. 

ftx 



(3* 



Equation (j3SJ) is rather difficult to study and we put one more restriction 
on /. Suppose that / = a(t) + b(ti) + c(t x ). Find the variables p, q, Yi(q), 
K_i(p) in terms of a, b, c 

p=-D- 1 (f tl ) = -b'(t), 

ft a'(t) 



Y 1 (q) = (d t + b'(t)d tx ) 



ft, d{t x ) ' 

a'(t) a"(t) b'{t)a'{t)c"{t x ) 



c(t x ) d{t x ) (c'(t x )) 2 



Y-i( P ) = (d t - ^-A x ){-b\t)) = -b"(t). 



c'(t 

Substitute these expressions into (T3"3"j) 



D x G(t,t x ) 



a"(t) _ V(t)a'(ty'(t x 
c^Q (C(t x )) 2 



where, 



G(t, t x ) 



1 ' C(t x ) + (c'(t x .)) 2 



(39) 



a 1 (t) 



b\t) 



(40) 



Evidently the left hand side of the equation fl39l) is of the form ^jrt x + ^-t xx 
and it contains the variable t xx while the right hand side does not contain it. 
This gives an additional constraint — 0. 
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Investigation of the equation (I3U|) is tediously long. Thus we will only 
give the answers. 

Theorem 3.1. If the equation ^ with particular choice of f(t,ti,t x ) = 
a(t) + b(ti) + c(t x ) has the operators Y\ and Y-\ such that the Lie algebra 
generated by these two operators is two-dimensional then f(t,ti,t x ) is of one 
of the forms: 

1) f(t,t u t x ) =7(t x + /3sinh(at + A)) +/3coshti 

2) f{t, h,t x ) = 7 (t c + (3e at ) + (3e at ^ + V , 

3) f(t,h,t x ) = c(t x ) + 7*1 + ^, 

4) f(t,t 1 ,t x )=j\n\t x \ + -ln(e t -e)+/3, 

7 

5) f(t, h, t x ) = jt x 2 + (3t x + at + V , 

where c(t x ) is an arbitrary function and a, (3, 7, A ; 77 are arbitrary constants. 

Additionally, if in the cases 1), 2), 3), the corresponding characteristic Lie 
algebras are also two-dimensional then the equations are of the forms: 
i) tlx = t x + (3 sinh(t) + (3 cosh(ti), 
H) h x =t x + e f + e' 1 , 
Hi) t\ x t x , 

and they have the n-integrals respectively: 

a) I = 4- cosh 2 1 — j3t x cosh t + y + (3t x sinh t — t xx + — n, 

b) i = *L + £-t 

U J 1 2 2 L XXl 

c) I t xx . 

Note that ii) has also x- integral as 

F = e* 1 "* + e 2 * 1 "* 2 "* + e* 1 "* 2 . (41) 

Thus, ti x = t x + e* + e' 1 is a discrete analog of the Liouville equation. The 
proof of the theorem 13.11 is given in the next section. 
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4 The Proof of the Theorem I3TT 



Here we consider the cases which satisfy the condition (139]) and prove the 
theorem 13. 11 

Case 1. a' = 0: 

The condition fl39l) turns out to be 

D ^)= h "- (42) 

This gives us h' = 7 so &(ii) = yti + (3. There is no condition on c(t x ). Hence 
/ becomes 

f(t,t 1 ,t x )=c(t x )+Tt 1 + P, (43) 
where (3 and 7 are arbitrary constants and c(t x ) is an arbitrary function oit x . 

Case 2. b' = 0: 

The condition fl39l) turns out to be 

fa'" a' -{a") 2 



a' J V (a') 2 

a" 



7 

Since only c depends on we have c = — which implies 

c(t x ) =7 In 1**1 +P- ( 44 ) 
The remaining equation gives a(t) as 

a{t) = -\n{-l+e xt+Xa + 0- (45) 
Thus, in a simplified form / is 

f(t,t 1 ,t x )=jln\t x \ + -ln(e t -e)+(3, (46) 

7 

where j3 and 7 are arbitrary constants. 
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Case 3. d" = 0: 

Note that we should have 

9G , . . _. 

ir (t,t x ) = (47) 

in the condition ( 1391) since there is no t x explicitly in the right hand side. If 
we expand (j4"71) we obtain 

c , ''c'[b\a') 2 + c\b') 2 a'}-(cy[b\a') 2 + 2c\b'ya'} + (c') 2 c''[a''b , -a'b''} = 0. (48) 

Under the case d" = we have c" = A and naturally d = At x + B where A, 
B are constants, we have 

A 2 t 2 [a"b' - a'b"} + 2At x [B(a"b' - a'b") - A(b') 2 a'} 

+ [-Ab'(a') 2 - 2AB{b'fa' + B 2 a"b' - B 2 a'b"} = 0. 

The coefficient of t x 2 should vanish, so a"b' — a'b" = 0. We should also have 

[B(a"b' - a'b") - A{b') 2 a'} = (49) 

from the coefficient of t x . Since a"b' — a'b" = 0, the equation becomes 

A(b') 2 a' = 0. (50) 

To provide this equation to be satisfied we have three choices. 

Subcase 3.i. A = 0: 

In this case c" = and so d = B where B is constant. The condition ( |39l 
becomes 



D x 

Hence 



a" + Bb" 
a' + Bb' 



Bbb" 



a" + Bb" 

— = a, 52 

a' + Bb' K J 

where a is constant. Thus we have two unknowns and two equations which 
are 

a' = ae at - Bb', (53) 
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and 

{a'f = B\b'f + p, (54) 

where a and p are constants. We take square of the first equation and 
subtract form the second one. We get 

(7 2 e 2at _ 2aBb'e at - p = (55) 

whose solution is 

b = ^ cosh(at + log + 5. (56) 
i?a yfp 

In a similar way we find a as 

a = sinh(at + log — — ) + r. (57) 
Ba Jp 

Here p ^ and 5, r are constants. Since c' = B, we find c^) = Bt x + (. 
Finally, after appropriate transformation we obtain 

f(t,t u t x ) = ~/(t x + /3smh(at + A)) +/3coshti +77, (58) 

where a, (3, rj, 7 and A are constants. 

If p = 0, then the equation (155]) gives us 

6 = — e at . (59) 

By using this equation we obtain a as 

< 60 > 

We have also c(t x ) = Bt x + £. Thus simplified form of / is 

f(t, h,Q = i(t x + (3e at ) + (3e at ' + 77, (61) 
where a, (3, 7 and rj are arbitrary constants. 
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Subcase 3.ii. a' = 0: 



In the Case 1. we have analyzed this case. Here, 

f(t,t u t x ) = a(t) + 6(ti) + c(t x ) 

= lt x 2 + f3t x + ah + r], 

where a, ft, 7 and r\ are arbitrary constants. 
Subcase 3.iii. b' = 0: 

Here c{t x ) = •yt x 2 + (3t x + A, but for b' = we have the following equation 

f a"'a'-(a") 2 ^ 
x \a>) ~ \ (a') 2 ) x 
a" 
c 7 ' 

which cannot be satisfied by this c(t x ). Hence a" = 0. So we obtain, 

f(t,h,t x ) = a(t)+b(h) + c(t x ) 
= lt x 2 + (3t x + at + r], 

where a, (3, 7 and 77 are arbitrary constants. 

The functions given in these cases may not give us n-integral / even though 
they satisfy the condition Let us see an example to this fact. This 

example belongs to Case 1. 

Example. 

ti x = f(t,t 1 ,t x )=t 1 + t x . (62) 
By using the definitions of the operators Yj, we obtain 

Yi = d t + d tx +d txx +d txxx + ... 
Y 2 = d tx + 2d txx + 3d txxx + ... 
Y 3 = d tx + 3d txx + Qd txxx + ... 
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and so on. As we see that the set Y 1 ,Y 2 ,Y 3 , ...,Y n is linearly independent 
which means that the algebra is infinite dimensional. By the theorem 11.11 
n- integral does not exist. 

Now we will focus on three problems which yield from the above consid- 
erations. We will study on n-integrals of these problems. 

Problem 1. The equation for the function / is 

tlx = fit, tx, t x ) = Xt x + \(3 sinht + (3 coshtx + 7. (63) 

Hence, 

t- Xx = g{t,t x ,t-i) = -(t x - A/3sinht_i - (3cosht - 7). (64) 
To obtain the operators Y\ we find the equations written below: 

= /Jsinhtx, 

= \t xx + \/3t x cosh t + (3ti x sinh 1 1 , 

= (3 cosh t x (Xt x + X(3 sinh t + (3 cosh t x + 7) + (3 2 sinh 2 1\ , 

= ^t xxx + \f3{t x ) 2 sinht + X(3t xx cosht + (3{ti x ) 2 coshti 

+ (3 svahti{\t xx + \f3t x cosh.t + (3ti x sinh ti), 

= 3f3 2 (t Xx ) sinh t x cosh t x + (3(t lx ) 2 sinh t x + (3 3 sinh 3 t x 

+ (3 cosh.ti{\t xx + \(3t x cosht + fit Xx sinhtx), 



fti 

fx 
fxt\ 
fxx 



fxxti 
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fxxx — Ai rax + \f3t xxx cosh t + 3\f3t x t xx sinh t + \(3t x cosh t 

+ 3(3 coshti(tix)[Xt xx + cosht + f3t\ x sinhti] + (3(ti x ) 3 sinhti 

+ (3 smhti{\t xxx + X(3t xx cosht + A/3t 2 sinht + (3(t\ x ) coshti 

+ [Atrrr + \(3t x cosht + f3t\ x sinhti]/? sinhti}, 

/xxxti = 3/5[At ra + A/9t x cosht + sinhti] [/3 sinhti cosh ti + ^ sinhti] 

+ 3/3ii x cosh ti[/? 2 sinh 2 ti + (3t lx cosht x ] + 3/3 2 (ti x ) 2 sinh 2 ti 

+ /3 cosh ti [At^ + A/St^ cosh t + A/3t 2 sinh t + /^t 2 ^ cosh t\ 

+ /3 sinh ti(Xt xx + \(3t x cosht + (3t\ x sinh ti)] + (3t\ x cosh ti 

+ (3 sinhti[/5 coshti(At ra + \(3t x cosht + (3t\ x sinhti) + /3 3 sinh 3 ti 

+ 3f3 2 ti x cosh ti sinh ti + f3t\ x sinh t x ] , 



Thus, 

D-\f tl ) = /3sinht, 
D-\f xtl ) = (3t x cosh t + /3 2 sinh 2 1, 

D~ l (f xxtl ) = 3(3 2 t x sinh t cosh t + /3t 2 sinht + /5i xa; cosh t + /3 3 sinh 3 1, 

19 



D'Hfxxxti) = 4/? 2 t 2 sinh 2 t + 6(3% cosh t smh 2 1 + 4 (3 2 t xx sinht cosh t 
+ 3(3 2 t 2 x cosh 2 1 + (3t 3 x cosh t + 3(3t x t xx sinh t 
+ (3txxx cosh t + /5 4 sinh t, 

Let us define (3 sinht = ip. Hence we can write the operator Y\ as 

+ {[{^ 2 + i) x ) x + ^OA 2 + + i>l(i> 2 + *Px)x + Hi> 2 + ^x)]}d txxxx 

+ .... (65) 
Now we will obtain the operator Y_i. From (|64|) we have, 
9t-i = -/3cosh*_i, 

&c = [txx - A/3t_ la; cosh t_i - /% x sinh t}/\, 

9xt-i = 1 cosh 2 t-x - (3t_ lx sinh 

= [^xxx — cosh £_i (£3.3. — A/3£_i x cosh t_i — sinh i) 

— /3t I:r sinht — (3t 2 x cosht — A/3(t_ia;) 2 sinhi_i]/A, 

9xxt-i — [3A/3 2 (t_i x ) sinh t_ x cosh t-\ — \(3{t^i x ) 2 cosht_x — A/? 3 cosh 3 t_i 

— {3sxnht^i(t xx — \j3t_ lx cosht_! — f3t x smht))/\, 
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Qxxx = {txxxx + -3fit-. lx sinh t_i (t xx - A/3t_i x cosh t_i - (3t x sinh t) 

— \f3(t^i x ) 3 cosht„i — (3t xxx sinht — 3(3t x t xx cosht — (3t x sinht 

— (3 cosh t-i [t xxx — (3 cosh t-i(t xx — \f3t-\ x cosh t_i — (3t x sinh t) 

— A / 3(t_ 1:E ) 2 sinht_ 1 - (3t xx sinh t - (3t 2 x cosht]}/ A, 

9xxxt-i = {(^x — ^Pt-i x cosht_! — /3t x sinht)(4/3 2 cosht_i shinty — 3(3t_ lx cosh£_i) 

— 3f3t-i x sinh t_i(2A/? 2 cosh 2 t_i — \f3t-i x sinht_i) 

— /3 sinh t-i [t xxx — [3 cosh t-i(t xx — \(3t-i x cosh t_i — (3t x sinh t) 

— \j3(t_i x ) 2 sinht_! — A/3(t_ia;) 3 sinht_i — (3t xx sinht — /3t 2 cosht] 
+ \j3 A cosh 4 1_! + 4X(3 2 t 2 _ lx cosh 2 t_i}/A, 



From these equations we get 
Digt^) = -13 cosht, 

^(^xt_i) = cosh 2 1 - /3t x sinht, 

D(g xx t_ 1 ) = 3(3 2 t x sinht cosht — /3t 2 cosht — /5t x x sinht — /3 3 cosh 3 1, 
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D(g xxxt _ 1 ) = 4:(3 2 t xx sinht cosh t — Spt x t xx cosh t — 6(3 3 t x cosh 2 t sinht 
+ Ap 2 t 2 x cosh 2 1 + 3p 2 t 2 x sinh 2 t - pt\ sinh t 
— Pt xxx sinh t + p A cosh 4 t, 



If we define —f3 cosht = ip we write the operator Y_i as 

y_i = 9i + (pd tx + (y? 2 + f x )d txx + [(^ 2 + y»x) a + V(V 2 + fx)]d txxx 

+ {[(v? 2 + <p x )x + v(v 2 + <Px)]x + vl(v 2 + (Px)x + viv 2 + Vx)]}d txxxx 

+ .... (66) 

Now let us only take the first three terms of Y\ and Y_i and find n-integral / 
such that both equations Y\l = and Y-\l = are satisfied. The first three 
terms of Y\ are 

Y 1 = d t + P sinh td tx + (P cosh t x + p 2 sinh 2 t)d txx . (67) 
We will use the method of characteristic, 

dt ^ dt x dt xx ^ dl (rq\ 

7 ~ p sinht ~ (p cosh t x + p 2 sinh 2 t) ~ T' 1 ' 

By the first and second terms we have 

P cosh t — t x = C\ . (69) 

By the first and third terms we have 

Pt x sinht -t xx = c 2 - (70) 

Hence from Y\l = we get the solution I = F(ci, c 2 ) where F is an arbitrary 
function. Now by using the transformation t — t, C\ — P cosht — t x , c 2 = 
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f3t x smht — t xx , we will write Y_\ in terms of t, C\ and c 2 . Note that 

d t = + ~lTTr + -7T-TT = d i + (3sinhtd ci + f3t x cosht<9 C2 , 

d t d- t d t d Cl d t d C2 

d~ t d d Cl d d C2 d . ~ 

d t* = 7T7T + 7T7T + T^TT - = + ^ smh ^ C2 , 

$ = + + — = -9 C . 

dt xx dt dt xx d Cl d txx d C2 

Under these transformations Y_i becomes 

F_i = <9 t ~ + P(smht + cosht)(9 cl — /3ci (sinh t + cosh £)<9 C2 . (71) 

Again by using the method of characteristic, 

dt dc\ dc 2 dl 

1 /3(sinht + cosht) —/3ci (sinh t + cosh t) 

we obtain 

/5(sinht + cosht) — ci = Ci (73) 

and 

| + c 2 = c 2 . (74) 
Thus from Y-\I = we obtain the solution 

c 2 

/ = G(/3(sinht + cosht) - Cl ,-± + c 2 ). (75) 

Z 

But we search for the common solution of Y X I = and Y_ X I = 0. Since the 
solution of Yil = only depends on C\ and c 2 we take the solution I = H{u) 

c 2 

where u = — + c 2 . Note that we may simply take 

Z 

c 2 $ 2 t 2 

7 = ^ + 02 = ^ cosh 2 i - 0t x cosh t + + f3t x sinh t - t xx . (76) 

Z Z z 

Remark 4.1. H^e should also check that Dl = I to be sure that I is n- 
integral. 

A 2 2 A 2 /? 2 l2 A 2 /3 2 7 2 • , 

£>/ = — t 2 + — — cosh 2 * — + — + \ 2 (3t x smh t 

z\ z\ z\ z\ 

+ \jt x + A7 sinh t — \t xx — \j3t x cosh t. 
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To make DI similar to I , we have to choose A = 1 and 7 = 0. In this case, 

Q 2 

Dl = l- P -. (77) 

So I is not n-integral, but it helps to find n-integral. Really, 

8 2 t 2 d 2 

F = y cosh 2 1 - (3t x cosh t + + (3t x sinh t - t xx + y n, (78) 

solves the equation DF = F. 

Problem 2. The equation for the function / is 

h x = f(t, h,t x ) = X(t x + e*) + e' 1 + c. (79) 

Hence 

t-i x = git, t x , = -(t x - e* + Ae*- 1 - c). (80) 
To obtain the operator Y\ we find the equations written below: 

fx ^t xx 4~ Xt x e 4~ t\ x e 1 , 

/** = ^txxx + A(t x )V + \t xx e l + (Xt xx + \t x e l + t lx e h )e h + (t lx ) 2 e h , 
fxxt! = e* 1 (Xt xx + At x e' + t lx e^ ) + St^e 2 ' 1 + e 3 ' 1 + (^) V 1 , 
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fxxx A£ xxx 4~ 2Xt x t xx e 4~ Xt x e 4~ Xt x t xx e 

+ 3t lx (Xt xx + Af x e* + ti^e* 1 ) + (t^Je* 1 

+ e'^At^ + At^e* + At 2 e' + [Xt xx + Xt x e l + t Xx e^\e tx + (t 2 Je 41 }, 

/xxxt! = (A* M + At.e* + ti :E e tl )(3ti a; e* 1 + 3e 2tl ) + e'^^e* 1 + 3e 2tl ) 

+ e'^At^ + Xt xx e f + At 2 e* + [Xt xx + Xt x e l + t^e^Je* 1 + (tL)e* 1 } 

+ e'HfAtxx + A* x e* + t^e'V 1 + 3e 2tl * lx + e 3 * 1 + (i^V 1 } + [t lx fe l \ 

Hence 

-D 1 (fxxt 1 ) 
D (fxxxti) 

Now we are ready to write the operator Yi: 

^ = ft + e'ft, + (e*f x + e 2t )d txx + {t xx e l + 3t x e 2t + e 3t + t 2 J)d txxx 

+ (3t x t**e* + 7t 2 e 2 ' + 4i xx e 2i + 6t x e 3t + t 3 e* + t xxx e l + e 4 *)^™, 

+ ... (81) 
If we define t x + e* = ip, Yi becomes 
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t x c ~\- c , 

= t xx 6 ~\- 3t x C ~\~ € -\- t x € , 

3t x t xx & -\- 7t x & -\- 4t xx c -\- 6t x c -\- t x & -\- t xxx & -\- c , 



Yi = 9 i +e t ^+eV^+e'(^ 2 +^)9 t _+e t [(^ 2 +^)x+V'(^ 2 +^)]^_ + ... . 

(82) 

Now we will obtain the operator Y^\. From flHUj) we have, 
9t-i = -e*" 1 , 

S'x = [txx — t x e t — At_ia;e* _i ]/A, 
9xt-! = -t-ixe*- 1 + e 2t -\ 

9xx \txxx t x € t xx € (txx t x d t—lx& )c ^(^— lx) ^ ]/^) 

ffaat-x = e*- 1 [-(t xx - tj - Xt-tee*- 1 ) + SXt^e*- 1 - Ae 2 *- 1 - A(t_ lx ) 2 ]/A, 

Qxxx [txxx 3t x t xx G t x C ~\- 3t—\ x C (txx t x d Xt—\ X G ) A(t_i a; ) 6 

- e^{t xxx - t xx e l - ty - (txx - t x e l - At-i^e'- 1 )^" 1 - A(t 2 _ 1:c )e*- 1 }]/A, 

9xxxt^ = [Xe'-H-t^t-t^- 1 -Be 2 *" 1 ) 

+ - ^e* - At_ 1 ,e t - 1 )(-3t_ 1:l .e*- 1 + 4c 2 *- 1 ) 

6 \t xxx t xx e t x c (txx t x c Xt_i x e )e A(t_^ a ,)e } 
+ e 2 *- 1 (^ - ^e* - At-i^-V 1 + Ae 4 ^ 1 - A^fe'^/A 
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From these equations we get 

D{g xt ^) = e 2t -t x e\ 
D(g X xt-i) = e t xx + 3~t x e e (t x ) e , 

D(yQxxxt-\) 3t x t xx e 4~ 7i x e 4~ 4t xx e 6t x e t x e t xxx e 4~ e , 



Thus the operator Y_i is 

y_! = d t - e% x - (e% - e 2t )d txx + {-eH xx + 3t x e 2t - e 3t - (tjV)^ 



■r.r.i 



+ (St^e* + 7t 2 e 2t + At xx e 2t - 6t x e 3t - t\e l - t xxx e l + e 4t )d txxxx 

+ .... (83) 
If we define t x — e* — (p,Y-i becomes 

y_i = <9 t -e t ^-eV<9 t;c:c -e*(v9 2 +^)9 t _-e t [(v9 2 +^) :E +^(^ 2 +^)]<9 t:c _+... . 

(84) 

Now we will only take the first three terms of Y\ and Y_i and find the integral 
I such that both Y 1 I = and Y_iJ = are satisfied. The first three terms 
of Yi are 

Y 1 =d t + e t d tx + (e 2t + t x e t )d txx . (85) 
The method of characteristic gives us 

dt dt T dt TT dl 



(86) 

1 e* e 2 ' + t x e* v ; 

By the first and second terms we have 

e* - t x = ci. (87) 
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By the second and third terms we have 

t x e* - t xx = c 2 . (88) 

Hence from Y\l = we get the solution / = F(ci, c 2 ) where F is an arbitrary 
function. Now by using the transformation t — t, C\ — e* — t x , c 2 = t x e l — t xx , 
we will write Y_ x in terms of t, C\ and c 2 . Note that 

4 " d t d; + d t d Cl + d t d C2 ~ 1 C1 * ° 2 ' 

di d d c , d d C2 d .- 
dt * = fTfT + JT + 7T 7T = + e'<9 C2 , 



<h,. <h d txx d Cl <),,., <),,, 



Under these transformations Y"_i becomes 

Y_ x = d~ t + 2e*d ci - 2c ie '<9 C2 . (89) 
Again by using the method of characteristic, 

dt dci dc 2 di 



1 2e* -2cie* 

we obtain 



and 



(90) 



2e* - ci = ci (91) 



+ c 2 = c 2 . (92) 



Hence from Y"_i/ = we obtain the solution 

/ = G(2e*- Cl ,| + c 2 ). (93) 

But we search for the common solution of Y±I = and Y-iI = 0. Since the 
solution of Yil = only depends on c\ and c 2 we take the solution I = H{u) 

c 2 

where u = — + c 2 . Note that we may simply take 
2 

2 2t f 2 

I = ^ + c 2 = - + ^-t xx . (94) 
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Remark 4.2. We should also check that DI = I. 

A 2 A 2 c 2 
DI = —t 2 x + —e 2t + - + \H x e l + \ct x + Ace* - \t xx - Xt x e\ 

If we choose A = 1 and c = 0, we see that I becomes 

I= Y + f-t xx (95) 
and it satisfies DI = I i.e. I is integral. 

Now we will find the x-integral of (1791) with the above conditions. Hence we 
try to find F which satisfies, 

D x F(h, h,t) = F t2 t 2x + F tl t lx + F t t x = 0. (96) 

Here t lx = t x + e* + e* 1 and t 2x = t lx + e h + e' 2 = t x + e* + 2e h + e' 2 . We 
insert these into (I9"B"1) . we get 

F t2 (t x + e* + 2e 41 + e* 2 ) + F h (t x + e* + e* 1 ) + F t t a = 0. (97) 

To satisfy this equation the coefficient of t x and the other terms should vanish 
separately Thus we have 

F t2 + F tl + F t = 0, 

F t2 (e l + 2c' 1 + e* 2 ) + F tl (e* + e* 1 ) = 0. 
So we get the operators, 

Xi = 9 <2 + S tl + 9 t , 

X 2 = (e* + 2c' 1 + e t2 )d t2 + (e* + e* 1 )^ = 0. 
The commutator of these operators is [Xl,X 2 ] = X 2 . 

Now for simplicity, we may consider the operator X 2 = e~ t X 2 . Explicitly, 

X 2 = (1 + e* 1 -*)^ + (1 + 2e 41 -' + e t2 - l )d t2 . (98) 
We will make change of variables, 

T = t, 
T\ = t X - t, 
T 2 = t 2 - t. 
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The operator X 2 becomes 

X 2 = (1 + e T1 )d T1 + (1 + 2e T1 + e T *)d T2 . 

Now we use the method of characteristics; 

dri dr 2 



1 + e Tl 1 + 2e T i + e 7 " 2 ' 



Equivalently we have 



-de 



-de 



1 + e Tl 1 + 2e Tl + e 7 " 2 
Let e -7 " 1 = u and e" 7 " 2 = v. Hence we get 

du dv 



l + u v(l + i) + l 



which has the solution 



u 1 - 



u + 1 



c — 



u 



If we put e n = u and e T2 = v in this solution we find c as 

(e ri + 1)- 



c = e 



1 + 



e 7 " 2 



_P — g*l _i _|_ g2*l— *2~* _|_ g*l _ *2 



(99) 



(100) 



flOl^ 



Hence F is 

Thus, iij; = t x + e* + e* 1 is a discrete analog of the Liouville equation. 



(102) 
(103) 

(104) 
(105) 



Problem 3. The equation for the function / is 

h x = f(t, t u t x ) = r(t x ) + Ph. 

Hence 

t-ix = g(t,tx,t-i) = r~ x {t x - f3t). 



(106) 
(107) 



30 



To obtain the operator Y\ we find the equations written below 

fa = P, 

f x = [r{t x )] x + p{r{t x )+pt{), 

fxti — P" \ 

fxx = [r(t x )] xx + P[r(t x % + P 2 (r(t x )+ ph), 

Clearly the operator Y\ is 

Y 1 = d t + Pd tx + P 2 d txx + P% xxx + P% xxxx + ... . (108) 

Since the function g = r~ l (t x — fit) does not depend on £_]_, the operator Y_i 
is 

Y-x = d t . (109) 

Remark 4.3. Let us take the first three terms ofY\. Clearly, the operators 
Y 1 and K_i commutes. From Y_\I = d t I = we see that I is independent of 
t. From Y\I = we obtain I as 

i = f(t xx -pt x ). (110) 

Simply, we may take I = t xx — pt x . Now we check if DI = I. 

DI = t xx r tx (t x ). (Ill) 

If P = and r tx (t x ) = 1 the equation DI = I is satisfied. In this case the 
equation U06\) takes the form t\ x = t x and I = t xx is its n-integral. 

Acknowledgments 

The authors thank Prof. M. Giirses for fruitful discussions. One of the 
authors (AP) thanks the Scientific and Technological Research Council of 
Turkey (TUBITAK) and the other (IH) thanks (TUBITAK), the Integrated 
PhD. Program (BDP) and grants RFBR # 05-01-00775 and RFBR # 06- 
01-92051-CE_a for partial financial support. 



31 



References 



[1] G. Darboux, Legons sur la theorie generale des surfaces et les applica- 
tions geometriques du calcul infinitesimal, T.2. Paris: Gautier-Villars 
(1915). 

[2] A. M. Grundland, P. Vassiliou, Riemann double waves, Darboux 
method and the Painleve property. Proc. Conf. Painleve transcen- 
dents, their Asymptotics and Physical Applications, Eds. D. Levi, P. 
Winternitz, NATO Adv. Sci. Inst. Ser. B Phys, 278, 163-174, (1992). 

[3] V. V. Sokolov, A. V. Zhiber, On the Darboux integrable hyperbolic 
equations, Phys. Lett. A, 208, no:4-6, 303-308 (1995). 

[4] I. M. Anderson, N. Kamran, The variational bicomplex for hyperbolic 
second-order scalar partial differential equations in the plane, Duke 
Math. J. , 87, no:2, 265-319 (1997). 

[5] A. B. Shabat, R. I. Yamilov, Exponential systems of type I and 
the Cartan matrices, (In Russian), Preprint, Bashkirian Branch of 
Academy of Science of the USSR, Ufa, (1981). 

[6] A. N. Leznov, V. G. Smirnov, A. B. Shabat, Group of inner sym- 
metries and integrability conditions for two-dimensional dynamical 
systems, Teoret. Mat. Fizika, 51, no:l, 10-21 (1982). 

[7] V. E. Adler, S. Ya. Startsev, On discrete analogues of the Liou- 
ville equation, Teoret. Mat. Fizika, 121, no:2, 271-284 (1999), (En- 
glish translation: Theoret. and Math. Physics, 121, no:2, 1484-1495, 
(1999)). 

[8] F. W. Nijhoff, H. W. Capel, The discrete Korteweg-de Vries equation, 
Acta Applicandae Mathematicae, 39, 133-158 (1995). 

[9] B. Grammaticos, G. Karra, V. Papageorgiou, A. Ramani, Integrabil- 
ity of discrete-time systems, Chaotic dynamics, (Patras,1991), NATO 
Adv. Sci. Inst. Ser. B Phys, 298, 75-90, Plenum, New York, (1992). 

[10] I. T. Habibullin, Characteristic algebras of fully discrete hyper- 
bolic type equations, Symmetry, Integrability and Geometry: Meth- 
ods and Applications, no:l, paper 023, 9 pages, (2005) // 
arxiv : |nlin.SI/0506027| 2005. 



32 



N. Kh. Ibragimov, A. B. Shabat, Evolution equations with nontrivial 
Lie-Backlund group, Funktsional. Anal, i Prilozhen, 14, no:l, 25-36 
(1980). 

A. V. Mikhailov, A. B. Shabat, R. I. Yamilov, A symmetry approach 
to the classification of nonlinear equations. Complete list of integrable 
systems, (In Russian), Uspekhi Mat. Nauk, 42, no:4, 3-53 (1987). 

R. I. Yamilov, D. Levi, Integrability conditions for $n$ and $t$ de- 
pendent dynamical lattice equations, J. Nonlinear Math. Phys. , 11, 
no:l, 75-101 (2004). 

R. I. Yamilov, On classification of discrete evolution equations, Us- 
pekhi Mat. Nauk, 38, no:6, 155-156 (1983). 

V. E. Adler, A. I. Bobenko, Yu. B. Suris, Classification of integrable 
equations on quad-graphs. The consistency approach, Communica- 
tions in Mathematical Physics, 233, no:3, 513-543 (2003). 

M. Giirses, A. Karasu, Variable coefficient third order KdV 
type of equations, Journal of Math. Phys., 36, 3485 (1995) // 
arxiv : |solv-int/9411004| 

M. Giirses, A. Karasu, Degenarate Svinolupov KdV Systems, Physics 
Letters A, 214, 21-26 (1996). 

M. Giirses, A. Karasu, Integrable KdV Systems: Recursion Oper- 
ators of Degree Four, Physics Letters A, 251, 247-249 (1999) // 



arxiv : |solv-int/9811013 



M. Giirses, A. Karasu, R. Turhan, Nonautonomous Svinolupov Jor- 
dan KdV Systems, Journal of Physics A: Mathematical and General, 
34, 5705-5711 (2001) // arxiv : |nlln~SI/0101031| 



33 



